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CN ■ Abstract 
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' The increasing focus on unquenched lattice simulations has revived interest also in gluonic 

Tij" ■ screening masses, whose inverses characterise the longest length scales at which thermal 

■ fluctuations are correlated in a hot non-Abelian plasma. We fill an apparent gap in the 

■ literature concerning the theoretical structure of one of the relevant screening masses (the 
one which equals twice the Debye mass at leading order) , by showing that the next-to- leading 
order correction to it is perturbative and small. This surprising result appears to explain 
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semi-quantitatively why this particular channel yields the smallest gluonic screening mass at 
temperatures around a few hundred MeV (it couples to the energy density and to the real 
part of the Polyakov loop), even though it is not among the smallest screening masses at 
asymptotically high temperatures. 
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1. Introduction 



Screening masses, or inverses of equal-time correlation lengths, are a fundamental character- 
istic of the long-range properties of a thermal system. Indeed, the quantum numbers and the 
degeneracy of the excitation with the lowest screening mass indicate what kind of an effective 
theory it is that determines the infrared sensitive thermodynamic properties of the system, 
such as finite volume effects |Tj. In QED, for instance, correlators of magnetic fields display 
a vanishing screening mass, while correlators of electric fields reveal a non- vanishing "Debye 
mass"; this then shows that at the longest length scales only magnetic fields are significant 
in an Abelian plasma, and finite-volume effects are powerlike. 

In non-Abelian gauge theories such as QCD, it turns out that the situation with the 
screening masses is a bit more complicated than in QED. In fact, even the definition of what 
is meant by screening masses requires some care: electric and magnetic fields, on which our 
Abelian intuition is based, are no longer gauge-invariant objects. Because of these subtleties 
it was only in the mid-1990's that fully satisfactory gauge-invariant and non-perturbative 
definitions were given to colour-electric, colour-magnetic, and certain more refined classes of 
screening masses in a non-Abelian plasma |2]. 

Following the conceptual clarification of the gauge-invariant definition of gluonic screening 
masses in QCD, systematic lattice measurements could also be carried out in all relevant 
channels. We would like to mention, in particular, quenched lattice measurements in four 
dimensions [H H] ; unquenched lattice measurements via a dimensionally reduced effective field 
theory in three dimensions [5]; and, most recently, unquenched lattice measurements directly 
in four dimensions [6]. Of course, a systematic analysis of the same observables can also be 
carried out in the AdS/CFT framework [7j. 

The purpose of the present paper is to consider the screening masses within the weak- 
coupling expansion. A number of them fall into the general class of observables whose leading- 
order value is fixed by the Debye scale; this class includes also many real-time observables 
of current interest, such as heavy-quark diffusion and jet quenching. It has been found in 
several such cases that the next-to-leading order correction is large for phenomenologically 
interesting values of the gauge coupling [8]. Our results will produce a "counter-example" to 
this empirical observation, showing that it is also possible to find observables in which the 
next-to-leading order correction is small. 

2. General framework 

In order to implement the resummations that are needed for defining the weak-coupling 
expansion at high temperatures, we choose to view the screening masses with the help of 
the dimensionally reduced effective field theory for hot QCD |9j, called EQCD [10]. This 
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approach is certainly sufficient for clarifying the theoretical structure of the various screening 
masses and, at least on the semi-quantitative level, also for numerical estimates. The effective 
Lagrangian has the form 

C^ = \'^ [Ffj\ + TV [A, Aof + mlTx [AI] + ... . (1) 

Here Fij = {i/ g^)[Di,Dj], Di = di - ig^Ai, At = A'^T", Aq = A'^T", and T" are hermitean 
generators of SU(3). In three dimensions the dimensionality of is GeV. A 2-loop derivation 
of m|, g'^ in terms of the parameters of four-dimensional QCD can be found in ref. pjj. 

Correlation lengths are defined from the exponential fall-off of two-point functions of local 
gauge-invariant operators. Without a loss of generality we assume the two-point functions 
to be measured in the xs-direction. The independent channels can be classified according to 
discrete symmetries defined in the two-dimensional transverse (xi, 3:2)— plane. A particularly 
important symmetry is often called R, and corresponds to the CT-symmetry of the original 
QCD; in terms of eq. ([1]), it sets ^0 ~^ —Aq. "Colour-electric operators" are defined to be 
odd under this symmetry, while "colour-magnetic operators" are even |2j. 

With this notion, examples of operators from which colour-electric screening masses can 
be determined are Tr [j4o-Fi2] and Tr[j4Q]. In four-dimensional QCD, these correspond to 
ImTr[PFi2] and ImTr[P], respectively, where P is the yet untraced Polyakov loop (we 
assume that the center symmetry is broken in the "trivial" direction, as is certainly the case 
in the unquenched theory). Note that these two channels do not couple to each other because 
of a different parity in the transverse plane. Typical operators from which colour-magnetic 
screening masses can be determined are Tr [Aq] and Tr[F^^2]) but any other gauge-invariant 
local singlet operator such as the energy density works as well. In four-dimensional QCD, 
Tr [Aq] corresponds to ReTr [P]. 

It is important to note that in principle the operators Tr [Aq] and Tr [-^^2] couple to each 
other [12^ I13j . In other words, if we measure a correlation matrix between these operators, 
then the matrix includes non-diagonal components. It is possible, however, to diagonalize 
the correlation matrix at large distances, i.e. to find two orthogonal eigenstates which display 
different screening masses (see, e.g., refs. |14]). It is these eigenvalues of the two-dimensional 
Hamiltonian that we refer to as M2 and M3 in the following. In practice, the coupling between 
Tr [Aq] and Tr [F^2\ is very weak, both parametrically [13] and numerically [5], so it appears 
to us that it should play no actual role in our analysis. 

Assuming that ^ 9e/'^i indeed the case at very high temperatures (in which 

limit rriE ~ gT{Nc/3 + A^f/6), g'^ ^ g^T, where 5'^/47r = Og is the strong gauge coupling, 
Nc is the number of colours, and A'f is the number of massless quark flavours), we can view 
Aq as a heavy field and write down the parametric forms of various screening masses within 



2 



the heavy-mass expansior0. In particular, the smallest screening mass in the colour-electric 
channel, coupling to Tr [^0-^12]) has a well-known logarithmic term at the next-to-leading 
order [15], and the general form [2] 
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Ml «mE + ^(ln^ + ci ) , (2) 



where ci ~ 6.9 for A'^c = 3 [16]. This expression works reasonably well down to low tem- 
peratures, overestimating the "exact" value within EQCD by a modest amount [5j. In the 
colour-magnetic channel, we can expect the mass coupling to Tr [^q] to have, in the heavy- 
mass limit, the form 



M2«2m, + ^(ln^ + C2) . (3) 
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Roughly, the correction here represents a three-dimensional bosonic analogue of the binding 
energy of a heavy quark-antiquark system, like J/ip. As far as we can see it is a non-trivial 
fact, following from the analysis in sec. [H that the coefficient of the logarithm in eq. ^ 
agrees with that in eq. The colour- magnetic screening mass which is the smallest at 
asymptotically high temperatures can, in contrast, be obtained from the theory from which 
Aq has been integrated out [13j; it couples dominantly to Tr [^12] ^^'^ has the form 
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M3«^XC3, (4) 



where C3 ^ 10.0 for TVc = 3 [TT] . 



3. Non-relativistic limit 

Our goal now is to estimate the coefficient C2 in eq. ([3]) which, to the best of our knowledge, 
remains unknown. The situation is quite similar to that in the case of the screening masses of 
fermionic bilinears, which we have studied previously in refs. [181 119j . The two adjoint scalar 
fields form a bound state, and a formal scale hierarchy exists between the heavy scalar mass, 
ttie; the relative momentum between the bound state constituents, p ~ {g^m^/Tr)^ ; and the 
binding energy, /S.E ~ ^I/tt, such that /rriE ~ Ai? (logarithms and numerical factors have 
been omitted; note that in terms of the four-dimensional coupling the scales are separated 
only by ~ (g'/vr)^/^). This scale hierarchy can be employed for constructing a set of effective 
field theories, perhaps ultimately a scalar analogue of PNRQCD [20]. As argued in ref. [18] . 
however, at the level of the correction of ©(^l/vr), the whole procedure simply amounts to 
solving the Schrodinger equation in a two-dimensional Coulombic potential for the s-wave 

^We stress that at this point the scale hierarchy ttie ^ ffi/''" serves only as a theoretical organizing principle 
for the computation; in practical estimates various group theory and numerical factors need to be added, and 
the phenomenological viability of the description can only be estimated a posteriori. 
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state; the only complication is that the heavy constituent "rest mass" entering the bound 
state problem needs to be fixed by a proper matching computation. 

To nevertheless give a somewhat more concrete indication of the effective theory setup, let 
us carry out a Wick rotation from the 3-dimensional Euclidean theory to a (2+l)-dimensional 
Minkowskian theory, and rename the xs-coordinate to be time, t. Let us, furthermore, write 
the time dependence of the quadratic part of the action in Fourier space, with uj denoting 
the frequency: 

Ao{uj,x)^ + ... . (5) 

If we concentrate on modes close to the on-shell points, |(^±mE| ~ sI/tt <C m^, and write lo = 
rriE+uj' or UJ = —m^—u)', then we observe that the dynamics of the forward-propagating mode 
ylQ(i-i;',x) = AQ{m^+uj' ,:x) and the backward-propagating mode (a;',x) = ^o(— — w', x) 
is determined by the non-relativistic Lagrangian 

In configuration space, the original field is related to the new effective fields by = 
e~*'"^*Ag -|- e™^*^Q . At leading order, then, the forward-propagating part of the composite 
operator Tr [AQ(t)] has the energy eigenvalue 2mE. 

When this argumentation is promoted to the quantum level, we expect the derivatives 
appearing in eq. ^ to get replaced by covariant derivatives, dtA'^ — > [L>t,ylQ]; the rapid 
oscillation frequency to get replaced by a matching coefficient, which we denote by 
Mrest! and the parameter in the denominator of the kinetic term in eq. ([6]) to be- 
come another matching coefficient, which we denote by Mkin- In the limit Mkin oo 
the propagators of the A'qS are replaced by Wilson lines in the adjoint representation: 
G(t,r) EE (^(f(t,r)A^'^(t,oX(0,rX(0,O)) = Tr {C/adj(t,r)C/Jdj(t,0)}, where C/adj(i,r) is 
a straight timelike adjoint Wilson line at spatial position r and we have for brevity omitted 
the (non-unique) spacelike connectors that make the point-split operators gauge- invariant. 
The evaluation of this expectation value leads to the concept of a static potential in the usual 
way: V{t) = \mit-,oo[idtG{t,r)\G^'^{t,Y). For the actual bound state problem Mkin ~ 
stays finite and the static potential takes the role of a matching coefficient. We do not need 
to know more about the effective theory setup in the following but remark that a formal 
discussion can be found in ref. 1221 . 
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Tr <^ Ao(-a;,x) 
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4. Determination of M2 



Proceeding now with the non-relativistic setup outlined above, we expect that in the heavy 
mass hmit the bound state mass can be written as 



M2 ^ 2Mrest + . 



(7) 



In dimensional regularization in d = 3 — 2e spatial dimensions, the next-to-leading order value 
of the matching coefficient Mrest reads [2] 



rest 



+ 1 



(8) 



where Ji is the scale parameter of the MS scheme. The binding energy can be solved from a 
two-dimensional Schrodinger equation; the potential appearing in it, obtained by integrating 
out the time (or X3) components of the gauge fields reads 



(2^) 



2-2e 



+ 27e 



In total, then, we are looking for the ground state solution to the problem 



2M, 



rest 



+ V{r) 



^0 = M-i"^ 



2^0 



where the non-kinetic terms combine to the finite expression 



2M,est + ^(r) = 2mE + 



2tt 



In(mEr) -|- 7e 



(9) 



(10) 



(11) 



In the kinetic term of eq. (jlOp . we already expanded the ("reduced" version of the) matching 
coefficient Mkin to leading order in g^/irm^, as is sufficient at the current level of accuracy. 

It is important to note that, unlike speculated in earlier works [21], no infrared divergences 
appear in eq. (jlip . The reason is that the logarithmic divergences originating from the "hard" 
momenta (g ~ m^] eq. ([8]); viewed from this side 1/e is an infrared divergence) and the "soft" 
momenta {q ~ 1/r; eq. ([9]); viewed from this side 1/e is an ultraviolet divergence) of the 
spatial gluons Ai cancel against each other in eq. (fTTI) . 

Carrying out suitable rescalings, eq. (jlOp can be solved up to one transcendental number. 
We thus obtain 

9lN, 



where 



M2 « 2mE + 



P = 



,0.60372466 - -Inp 
27r V 2 



gjNc 

27rmE 27ryi/2 
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(13) 



5 




Figure 1: The parameter m^/g^ — y^^"^ from eq. (fT4|) for Nc = 3. 

At next-to-leading order in massless QCD the ratio y = m'^/g^ is renormalization group 
invariant [23], and can be written compactly as 



y 



{2N, + N{)inN, - 2N{) 
1447r2 



' A-kT 4iVfln2-iVc 5N^ + N? + 9Ni/2Nc 

In 7e + -TTTT + 



A- 



MS 



lliVc-2iVf {2Nc + Nf){llNc - 2N{) 
The corresponding m^/g^ = y^l"^ is plotted in fig. [1] for = 3. 



(14) 



5. Summary and conclusions 

Comparing eq. (fT2]) with eq. (I3|), we obtain 



C2 1.9467141 



(15) 



for Nf^ = 3. Given that > 0.5g'^ (cf. fig.[I|), the latter term in eq. is always subdominant. 
This is in stark contrast to eq. ([2]), in which the latter term, containing the non-perturbative 
coefficient ci ~ 6.9, dominates in the whole temperature range of phenomenological interest. 
Note that because C2 <C C3 « 10.0, M2 is in general also below M3 (cf. eq. (jH)) in the 
temperature range of fig. [TJ All three masses are plotted in fig. ^ both in units of g^ and in 
units of T. 
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Figure 2: Left: The masses Mi,M2,M^ (cf. eqs. ©-([I])), in units of g^, for A^c = 3, after insertion 
of m-E/ from eq. (|14p . For better visibility, the four curves for M3, which are degenerate, have been 
sHghtly displaced from each other. Right: The masses M2,M3, which belong to the same quantum 
number channel, after the insertion of the 2-loop value for g^/T from ref. pT| . 

To summarize, it appears understandable that M2 represents the smallest screening mass 
at realistic temperatures, because of the small perturbative coefficient C2 ~ 1.9 in its next-to- 
leading order term, even though in the extreme limit ^ 9^/'^ it eventually overtakes both 
Ml and M3, because a higher multiple of rng's appears in the leading term. This observation, 
together with the explicit results in fig. [51 constitute the main points of this note. 

We would like to stress, finally, that although our result for M2 is not meant to be quan- 
titatively accurate at low temperatures, it nevertheless reveals an interesting pattern. For 
example, for A'^f = 0, fig. [2] suggests M2/T ~ 3. ..4 in the phenomenologically interesting tem- 
perature range, while lattice measurements indicate values M2/T ~ 2.5. ..3 Od], i.e. in the 
same ballpark but deviating downwards on the quantitative level. It appears, though, that 
this difference could at least partly be understood through higher order corrections within 
the EQCD effective theory defined by eq. ([T]): for T/A,^ 2, the non-perturbative lattice 
measurements in ref. [5] yielded M2I g\ ~ 1.0, 1.3, 1.5, 1.6 for A^f = 0, 2, 3, 4, respectively, and 
Mxlg\ K, 1.7,2.0,2.1,2.1 for the same cases; these values lie consistently somewhat below the 
perturbative estimates in the left panel of fig. [21 resulting in a better accord with 4d lattice 
data. Moreover, a similar overshooting of the C((7g)-corrected screening masses has been 
found for mesonic observables |18t I19j . So, it might be the general case that higher-order 
corrections, mostly from within three-dimensional EQCD dynamics, sum up to a negative 
correction to the next-to- leading order expression for screening masses. This would imply 
that screening masses are different in character from some dynamical quantities like the heavy 
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quark momentum diffusion coefficient, in which case higher order corrections appear to add 
up on top of the aheady large next-to- leading order correction [24] . 



6. An open issue 

We end by briefly pointing out an open problem to which we have no solution. Consider 
the screening mass extracted in four dimensions from the imaginary part of the Polyakov 
loop; in EQCD this corresponds to M(Tr[74Q]). In weak coupling, M(Tr[AQ]) = Suie + 
^^^(c4ln^ + C5). Even though this is heavier than Mi, particularly for A'^f > it 
can be easily measured on the lattice [6j, since the corresponding operators have different 
quantum numbers. Therefore, it would be nice to know the coefficients 04,05. Though 
this is certainly a well-defined problem (for Nc > 2), we have no clear idea about how it 
could be solved in a systematic way. (Probably the system can still be described by a non- 
relativistic many-body Schrodinger equation with a certain three-body potential in it, but in 
the absence of an effective theory framework or an explicit power-counting argument, it is 
difficult to know for sure how to proceed without ambiguities.) If the pattern found in this 
note continues, however, we might expect C4, C5 to be coefficients at most of order unity, such 
that the leading-order term would dominate even more than in M2. 
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